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Abstract. The observed hard X-ray (HXR) flux spectrum 7(e) from solar flares is a combination of primary 
bremsstrahlung photons 7p(e) with a spectrally modified component from photospheric Compton backscatter of 
downward primary emission. The latter can be significant, distorting or hiding the true features of the primary 
spectrum which are key diagnostics for acceleration and propagation of high energy electrons and of their en- 
ergy budget. For the first time in solar physics, we use a Green's function approach to the backscatter spectral 
deconvolution problem, constructing a Green's matrix including photoelectric absorption. This approach allows 
, spectrum-independent extraction of the primary spectrum for several HXR flares observed by the Ramaty High 

Energy Solar Spectroscopic Imager (RHESSI). We show that the observed and primary spectra differ very sub- 
stantially for flares with hard spectra close to the disk centre. We show in particular that the energy dependent 
photon spectral index 7(e) = — d log 7/dlog e is very different for 7p(e) and for 7(e) and that inferred mean source 
electron spectra F(E) differ greatly. Even for a forward fitting of a parametric F(E) to the data, a clear low-energy 
cutoff required to fit 7(e) essentially disappears when the fit is to 7p(e) - i.e. when albedo correction is included. 
The self-consistent correction for backscattered photons is thus shown to be crucial in determining the energy 
^ ' spectra of flare accelerated electrons, and hence their total number and energy. 

> 

X 

C3 , 1. Introduction observed spectrum may have been substantially affected 

by backscatter. This effect is well-known in solar physics 

Bremsstrahlung HXR spectra of solar flares give us direct ( and more generally in X-ray astronomy - Magdziarz and 

information about the distribution of accelerated high en- zdziarski 1995) and there have been several discussions of 

ergy electrons propagating in dense regions of the solar its i nnuen ce on observed X-ray spectra (Tomblin, 1969; 

atmosphere. Though quite dense, the atmosphere above Santangelo et al, 1973, Bai and Ramaty, 1978) , and on 

the X-ray bremsstrahlung producing region is optically electron spectra inferred from them (Johns and Lin, 1992; 

thin, but the very dense lower photospheric layers are not Alexander and Brown, 2002), now that high resolution 

and photons emitted downwards are Compton backseat- data are available. Contamination of the observed X-ray 

tered by atomic electrons there with an energy depen- spe ctrum by reflected photons leads to flattening of the 

dent efficiency. Photons of energy above ~ 100 keV pen- spe ctrum and underestimation of the spectral index (Bai 

etrate so deeply that they are lost to the observer, while and Rama ty, 1978). This can lead to substantial under- 

at below - 10 keV they are mostly photoelectrical^ ab- estimation of the total extrapolated electron energy in a 

sorbed by their first scattering. The resulting spectral re- nare . Spectral features connected with a low energy cut- 

flectivity thus has a broad hump in the range 10-100 keV off in the electron distribution, of crucial diagnostic inter- 

with a maximum around 30-40 keV. At some energies est for determining the acceleration mechanism and the 

and view angles the reflectivity approaches 100% so the total fast electron energy content, are hidden or substan- 

tially changed by the backscattered component (Zhang 

Send offprint requests to: Eduard P. Kontar and Huang, 2004). Moreover, an apparent low-energy cut- 
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off in albedo-uncorrected electron spectrum may not be 
required by data from the "true" primary (albedo cor- 
rected) spectrum. 

Comptonisation of prescribed primary HXRs can be 
straightforwardly studied using Monte-Carlo simulations 
(e.g. Bai and Ramaty, 1978). This technique is ideal for 
obtaining the reflected, and hence the total observed, pho- 
ton spectrum if one knows or assumes the form of the 
primary X-ray spectrum (Alexander and Brown, 2002). 
However, the primary spectrum is unknown and unlikely 
to be an exact power-law, as is usually assumed. Even 
with a power-law distribution of primary electrons a num- 
ber of physical processes exist which will cause the re- 
sultant photon distribution to deviate from pure power- 
law form e.g. effects of nonuniform ionization (Brown, 
1973; Kontar et al., 2003), beam return current losses 
(Zharkova and Gordovskyy, 2005) and others. Results ob- 
tained from Monte Carlo simulation for power-law pri- 
mary photon spectra cannot be straightforwardly applied 
to such situations. More sophisticated parametric mod- 
elling of the primary photon spectrum, with more theo- 
retical uncertainties and free parameters, would become 
necessary. Therefore, an approach independent of the pri- 
mary spectrum is required. 

Understanding and modelling backscatter has become 
even more important with the advent of high quality X-ray 
spectra from the Ramaty High Energy Solar Spectroscopic 
Imager (RHESSI) with spectral resolution as high as ~ 1 
keV, in combination with uncertainties as low as a few 
percent (for strong flares) allowing detailed spectral anal- 
ysis over a broad X-ray range (Lin et al, 2002). A correct 
treatment of Compton backscatter by a cold medium has 
also been important for understanding spectra from cos- 
mic X-ray sources (reviewed e.g. by Svensson, 1996). This 
led to construction (Magdziarz & Zdziarski, 1995) of an 
analytic (piecewise continuous), angle dependent, Green's 
function G(/i,e ,e) for the probability that a photon, of 
initial energy e Q and drawn from a distribution isotropic 
in the downward hemisphere, incident on a plane-parallel 
atmosphere, is re-emitted in direction (/z = cos 9) with 
energy e. This was obtained by fitting the results of a large 
set of Monte Carlo calculations, and its form is guided, 
particularly in various limits, by analytical results such as 
those of Illarionov et al, (1979), Lightman et al. (1981), 
White et al. (1988). Representing as it does the results of 
Monte Carlo calculations, no assumptions of single scat- 
tering or on the form of primary spectrum are involved. 
The only major assumptions in the Green's function for- 
mulation are those already mentioned: plane geometry of 
the cold matter - this is well satisfied for the photospheric 
conditions and downward isotropy of the primary photon 
distribution. 

Combined with a treatment of absorption (Morrison 
and McCammon, 1983) with standard solar photosphere 
element abundances (Anders and Grevesse, 1989), the 
Green's function formulation of Compton backscatter pro- 
vides us with a powerful tool for interpreting flare hard X- 
ray spectra. Firstly, it allows us to calculate easily, without 



time consuming Monte Carlo simulations, the total ob- 
servable photon spectrum for any specified primary spec- 
trum, not just the power-law cases treated by Bai and 
Ramaty (1978). More importantly, because the Green's 
function may be used to express the observed spectrum 
in terms of any primary spectrum, we may use regular- 
ized inversion techniques (Craig and Brown 1986, Kontar 
et al, 2004, and references therein) to deconvolve the pri- 
mary photon spectrum from observations. Features of this 
primary spectrum arise only from features of the HXR 
source electron distribution, not of the energy-dependence 
of backscattering, so it may be safely inverted (Kontar et 
al, 2004) to yield a meaningful source-averaged electron 
distribution in the sense of (Brown et al., 2003). Although 
the Green's function makes such model-independent in- 
version possible, this step is carried out here for the first 
time. 

Section 2, 3 describe the Green's function formula- 
tion and the incorporation of absorption respectively. 
Section 4 describes how to recover the primary spectrum 
from data, and Section 5 applies this to real RHESSI data. 
In Section 6 we discuss applications and implications of 
our work. 



2. Backscattering of X-rays 

As discussed in Section Q downward-emitted photons will 
be either absorbed or scattered, and some scattered to- 
ward the observer, adding to the total flux detected. The 
backscattered photons will be emitted from an extended 
region of the photosphere, out to the horizon distance but 
with greatest intensity in an area of extent comparable 
to the altitude of the primary source (Brown, McClymont 
and van Beek, 1975; Bai and Ramaty, 1978; Schmahl and 
Hurford, 2004). The Green's function implemented here 
ignores spatial structure, however, simply summing over 
all photons that emerge in the correct direction. This 
makes it ideal for dealing with spatially unresolved spec- 
tra. 

While the Green's function (Magdziarz & Zdziarski, 
1995) was obtained for use in studying cosmic X-ray 
sources, it applies to any plane-stratified slab of 'cold' 
material (i.e. the scattering electrons are non-relativistic) . 
Scattering takes place only on electrons, whether free or 
atomic. To account for elements heavier than hydrogen the 
Compton cross-section is multiplied by an effective mean 
atomic number z=1.2. The detailed density structure of 
the medium is irrelevant (Tomblin, 1972). 

Absorption, on the other hand, does depend strongly 
on chemical composition, and the best estimate of photo- 
spheric abundances should be included when adding ab- 
sorption to the Green's function treatment. We use stan- 
dard abundances given by Anders and Grevesse (1989). 
The heavy elements Fe/Ni play the most important role in 
the range from 10 keV up to 30 keV, while lighter elements 
contribute below 10 keV (Morrison and McCammon, 
1983). 
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3. Green's function 

The definition of the Green's function allows us to write 
the secondary, backscattered spectrum Is for any isotropic 
primary spectrum Ip (both have units of photons keV -1 
s" 1 cm -2 ): 

e moI 

Is(e)= J I P G(fi,e,e )de (1) 

e 

where G(fx, e, e ) is the angle-dependent Green's function 
and fj, = cos(#), where 8 is the heliocentric position angle 
of the X-ray source. The observed spectrum, at direction 
(U, is then the sum of Is and Ip. The importance of not 
averaging over viewing angle may be seen from Bai and 
Ramaty (1978). The Green's function itself is shown in 
Figure la - for details see Lightman and White (1988), 
Magdziarz & Zdziarski (1995), Ross and Fabian, (1993) 
and references therein. 

The total observed spectrum /(e) is given by the sum 
of Ip and Is, i.e. using Q 

I(e)=I P (e)+ J IpG(n,e,e )de (2) 

With measured 1(e) the primary spectrum Ip, which tells 
us directly about the source electrons, may be obtained 
via solution of the integral equation J5| • This is a Volterra 
Equation of the second kind (with Ip appearing outside 
as well as inside the integral, solution of which is more 
stable against noise in I than those of the first kind (e.g. 
Craig and Brown, 1986)). In addition, the kernel (Green's 
function) is peaked so that the kernel function/matrix is 
not too far from diagonal. Thus regularization to achieve 
an accurate stable solution for Ip given / is less trouble- 
some than in other inverse problems such as using Ip(e) 
to find the source electron spectrum following Kontar et 
al (2005). 

In practice measurements yield discrete quantities and 
the integral equation J5J) is presented in matrix form 

I(e i ) = I P (e i )+GMlp(ej) (3) 

where we have used the summation convention for re- 
peated indices, and introduced what we may call the 
Green's matrix 

Gij(fi) = J G(fi,e ,£i)de (4) 

Note, that the integration in Eq. Q is best performed via 
a change of variable to the wavelength domain y = mc 2 /e, 
j/o = mc 2 /eq due to sharp features in the Green's function 
(Figure Special care should be taken when integrating 
near Ay = y - y Q = 2. 

The shape of the Green's function depends on the 
maximum energy of a primary photon. For primary pho- 
tons with low energies eo < 30 keV, the Green's func- 
tion has a rather simple structure close to a Dirac's delta 



function (Figure ^) showing that backscattering is domi- 
nated by the first scattering, while the contribution from 
higher orders of scattering is small. In the range be- 
low ~ 12 keV photoelectric absorption dominates scat- 
tering and Compton scattering can be treated as diffuse 
monochromatic radiation transfer. The Green's function 
then can be replaced by C?(/i, e, eo) ~ F(/i, A)<5(e — eo) 
(Chandrasekhar, 1960, where the form of F(fi, A) is given), 
A(e) = o- c (e)/(<7c(e) + <Xyi(e)), o c {e), <r A {e) are the to- 
tal Compton and absorption cross-sections respectively. 
A gives the single-scattering albedo modified by absorp- 
tion and completely controls the influence of absorption. 
Application to the solar case simply requires evaluation 
of A with solar photospheric abundances. In this energy 
range below 10 keV the total Compton scattering cross- 
section can be approximated by the Thompson cross- 
section 87rr 2 /3. The absorption cross-section is based on 
the element abundances from Anders and Grevesse (1989). 
We use the coefficients of a piecewise polynomial fit to 
the numerical results given by Morrison and McCammon 
(1983) in the range below 10 keV and extrapolate them 
as a a (e) ~ e~ 3 above 10 keV. Due to rapid decrease of 
the absorption cross-section a a ~ &C around 12 keV. The 
generated Green's matrixes Q have uncertainty around 
6% (due to the approximations used) and 1 keV energy 
resolution. 

The absorption is photo-electric in nature and below 
10 keV has a complex structure with multiple absorption 
edges due to elements heavier than hydrogen. In conse- 
quence it is also strongly influenced by elemental abun- 
dances. Figure |21 shows the two absorption edges of Fe at 
7.1 keV and Ni at 8.3 keV in the reflected component. The 
spectrum shows features in agreement with the study of 
the reflection spectra in intergalactic nuclei by Ross and 
Fabian (1993). Here we ignored the fluorescence of the 
lines which is beyond the scope of our paper. 

Note that Bai and Ramaty (1978) and Zhang and 
Huang (2004) assumed a different absorption approxima- 
tion from us based on earlier photoelectric absorption 
cross-sections by Fireman (1974) and photospheric abun- 
dances by Withbroe (1971). However, we should note that 
the difference is rather modest accounting for around 30% 
of reflectivity at low energies and becoming smaller at 
higher energies, where absorption is less important (Figure 
|5J). The reflectivity, and thus the contribution of Compton 
back-scattered photons into observed spectra is spectrum- 
dependent (FigureOJ. Therefore, the albedo spectrum de- 
pends on the shape of the primary spectrum (Figure , 
- previous studies (Bai and Ramaty 1978, Alexander and 
Brown, 2002) considered only the results for prescribed 
power-law primary spectra. This difference is important 
when one tries to infer a primary spectrum from an ob- 
served one. 

We verified that we can reproduce the forward re- 
sults of Bai and Ramaty (1978) by numerically evaluat- 
ing Eq. Q for the particular case Ip ~ e~ 7 and the pa- 
rameters described in Figure[5] Consistency with previous 
treatments being established in this way, we now turn to 



4 



Eduard P. Kontar et al.: Compton backscattered X-rays from solar flares 



3 



1 00.00 f 



10.00 : 



1 .00 : 



0.10 r 



0.01 I I £_ 



1 00 



Energy, keV 



Fig. 1. Green's functions G(/x, e,eo) including Compton 
scattering and photoelectric absorption for three primary 
photon energies e = 20,80,300 keV and fx = 0.7 (9 ~ 
45°). 
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Fig. 3. Spectrum dependency of the reflected flux. Upper 
panel: Primary (dashed line), reflected (solid line) and 
total (dotted line) photon spectra calculated for a mod- 
elled flare spectrum Ip(e) ~ e -3 for e > 45 keV, Ip{e) ~ 
e~ 2 + exp(— e/1.7) for e < 45 keV. The lower panel shows 
the reflectivity for a X-ray source at heliocentric angle 
9 = 45°, defined as the ratio of reflected to primary fluxes 
R(e,9 = 45°) = I R (e)/I P (e) (solid line). The reflectivity 
assuming power-law for a primary source with spectral 
index 3 (dotted line) and 2 (dashed line). 
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Fig. 2. Upper panel: Primary (solid line), reflected (dotted 
line) and total (dashed line) photon spectra calculated 
assuming a primary spectrum Ip(e) ~ e~ 3 , and using the 
Green's function for a X-ray source at heliocentric angle 
9 = 45°. The lower panel shows the reflectivity, defined 
as the ratio of reflected to primary fluxes R(e, 9 = 45°) = 
Iji(e)/Ip(e). The reflectivity taken from Bai and Ramaty 
(1978) is shown with diamonds. Two absorption edges of 
Fe at 7.1 keV and Ni at 8.3 keV are seen in the reflected 
component. 



the inverse problem - deduction of the primary spectrum 
from observations of the total spectrum. 



4. From observed to primary X-ray spectrum 

In the absence of data noise we could immediately invert 
Eq.|21to obtain Ip: 



I P (e j ) = A ij I(e i ), A=(l + G(/i))- 



(5) 



The same result as Equation J5J can be obtained solving 
the original system of linear equations Q using, for exam- 
ple, singular value decomposition. In reality, we have pho- 
ton spectra contaminated by noise. If the resulting matrix 
has a large condition number (ratio of largest to small- 
est singular values), then a regularized approach should 
be adopted (e.g. Craig and Brown 1986, Piana et al 2002, 
Kontar et al, 2004). In practice, we can incorporate the 
effects connected with solar albedo as a part of instrument 
response. 

4.1. SPEX and albedo correction 

For a given instrument such as RHESSI, the count flux 
C(ei) is the linear composition of instrument response ma- 
trix Rjj and the photon flux coming from the Sun 
(Schwartz et al, 2002) 

C{ei) = Ry lie,) 



(6) 
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However, since we are interested in the flare primary spec- 
tra rather than the observed flux, we can combine equa- 
tions © and © to get 

R' 



17-Sep-2002 05:50:47.99 



(7) 



where R'y = R,-fc(ljy + aG^). Here we introduced a co- 
efficient a that accounts for anisotropy of the source in the 
Eddington approximation - that is, distinct but constant 
specific intensities in the downward and in the upward 
hemispheres, in the ratio a. 

This approach has already been incorporated into the 
SPEX software (Schwartz et al, 2002) including the ob- 
ject oriented version (OSPEX). The software is publicly 
available as a part of Solar Soft Ware (SSW). 

5. RHESSI solar flare HXR spectra 

5.1. Photon spectral correction 

The spectral effect of Compton back-scattering on HXRs 
is most pronounced for harder (flatter) spectra such as 
in the flare of September 17, 2002 around 05 : 50UT 
(see Fig. fJJ which we study here though quite a num- 
ber of such flares have been observed: August 20, 2002 
(Kasparova et al, 2005), April 25, 2002 (Kontar and 
Brown, 2005). The location of our flare, which uniquely 
determines the Green's function, was in the western hemi- 
sphere at (560, —300) arcseconds corresponding to a helio- 
centric angle of 9 ~ 41°. 

We consider the spectrum around flare maximum 
(5:50:48-5:51:36 UT ) when we have enough counts over 
the range 3 keV up to > 100 keV. We used 7 out of 9 
front segments, excluding detectors 2 and 7 due to their 
low energy resolution at the time of observation (Smith et 
al, 2002). The time interval is shown in Figure 4. 

After background subtraction, spectral results are 
shown in Fig. [5] for both the total observed spectrum, 
and the primary spectrum Ip determined via the method 
of the previous section and using a = 1 (isotropic emis- 
sion). We note how removal of backscattered photons, as 
expected, reduces the primary flux and steepens the spec- 
trum in the 20 - 50 keV photon energy range and slightly 
flattens the spectrum at energies above the albedo maxi- 
mum. 

5.2. Effect of albedo on photon spectral index 

The albedo changes the spectral indices of photon spectra, 
both of power-law fits and of the local energy dependent 
spectral index defined by 

dl(e) din 1(e) 



7 (e) ee - 



(8) 



/(e) de dine 

Direct numerical calculation of any derivative gives poor 
results in the presence of data noise, even for data as good 
as that from RHESSI. Differentiation is a specific case of 
solving an integral equation and to obtain meaningful re- 
sults in the presence of noise needs some sort of stabiliza- 
tion such as Tikhonov regularization, which we use here 
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Fig. 4. Temporal variation in four energy channels (4 sec- 
onds cadence) of the count rates summed over seven front 
RHESSI detectors for the September 17, 2002 solar flare. 
The vertical lines show the accumulation interval selected 
5:50:48-5:51:36 UT for further spectral analysis. 
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Fig. 5. Spectrum of the September 17, 2002 solar flare in 
the time interval 5:50:48-5:51:36 UT. The solid line shows 
the observed spectrum; the dash line is the primary spec- 
trum (corrected for Compton back-scattering). 

following Kontar and MacKinnon (2005). We limited our 
analysis to the most interesting range between 10 and 40 
keV, where the albedo effect is most pronounced. At the 
energies above 40 keV the count rate is small and thus 
uncertainties are large, whereas at the energies below 10 
keV Fe/Ni lines affect the bremsstrahlung spectra. 

Figure jnj) shows the energy dependent spectral index 
for our flare spectrum. Both spectra (with and without 
albedo correction) increase toward the ends of the energy 
interval. At the low energies this is probably due to a 
thermal component while at higher energies it is likely 
connected with softening of a nonthermal electron spec- 
trum. The spectral index of the observed data shows a 
clear minimum around 15 keV, while the spectrum cor- 
rected for albedo has a rather extended minimum in the 
broad range 15 — 30 keV. The albedo correction increases 
the minimum value of 7 from 2.0 to 2.6. The opposite ef- 
fect can be seen at higher energies: the spectral index at 40 
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Fig. 6. Energy dependent photon spectral index 7(e) of 
the September 17, 2002 solar flare. The solid line shows the 
spectral index of the observed spectrum; the dash line is 
the spectral index for the primary spectrum (corrected for 
Compton back-scattering) . The confidence intervals repre- 
sent the range of solutions found by allowing the incident 
photon spectrum to range randomly within the estimated 
(instrument + shot noise) errors. 

keV decreases from 3.5 to 3.2, albedo correction making 
the primary spectrum harder. Qualitatively similar results 
have been obtained for forward modelled spectra (Bai and 
Ramaty, 1978; Zhang and Huang, 2004). 

5.3. Effect of albedo on inferred source averaged 
electron spectra 

The primary photon flux is directly connected with the 
mean source electron (so-called thin target) spectrum 
(Brown, 1971; Brown et al, 2003) 

^( £ ) = i^2^/ F{E)Q{e,E)dE, (9) 

where Q(e, E) is the bremsstrahlung cross-section differ- 
ential in e (Haug, 1997) and the density weighted mean 
radiation source electron flux spectrum F(E) (electrons 
cm~ 2 s _1 keV -1 ) is defined by 

j?(£) = _L I F(E,r)n(r)dV. (10) 

with F(E, r) and n(r) the local electron spectrum and 
source proton density at position r in radiating volume V 
with mean target proton density n — V" 1 J n(r) dV. 

Solution of Equation @ for F(E) is a somewhat un- 
stable inverse problem (Craig and Brown 1986, Piana et 
al 2003) and so results are strongly affected when the ob- 
served /(e) is used instead of the albedo corrected Ip(e) 
required for correct solution of ©. F(E) is important in 
the study of acceleration mechanisms and of the flare elec- 
tron energy budget. Of particular interest are data such 
as from the flare discussed above which exhibit very flat 
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Fig. 7. Mean electron flux spectra of the September, 2002 
solar flare. The solid/dash lines show the spectrum with- 
out and with albedo correction. 

regions which may correspond to F(E) with a low energy 
cut-off in the mean electron spectrum which would have 
major implications for electron acceleration and propa- 
gation, and for the energy budget (Kontar and Brown, 
2005). 

We used standard SPEX software (Schwartz et al., 
2002) to fit to the observed X-ray counts an approxi- 
mate parametric F(E) comprising an isothermal plus bro- 
ken power-law form with low energy cut off . The mini- 
mum x 2 fit of the observed spectra gives emission measure 
EM = 8.0 x 10 46 cm" 3 , temperature kT = 1.36 keV, spec- 
tral indexes Si ow = 1.11, Shigh = 3.01 with break energy 
Eb — 54 keV and a low energy cut-off at 18 keV very 
clearly visible above the isothermal part (Figure EJ). 

Next we carried out the same model fit but using the 
instrument response matrix corrected for albedo as de- 
scribed above - i.e. we repeat the process for Ip instead 
of /. The result shows the same tendency as we saw in 
the analysis of energy dependent photon spectral index. 
The primary photon spectrum is softer at low energies and 
harder at high energies. Now we find a slightly bigger emis- 
sion measure EM = 9.2 x 10 46 cm~ 3 and slightly smaller 
temperature kT — 1.29 keV. The spectral index of ener- 
getic electrons becomes softer at low energies Si ow = 1-28 
and harder at high energies Shigh = 2.73 with the break 
energy at higher energies Eb — 59keV. Most importantly, 
in contrast with the fit to the observed spectrum, the pri- 
mary spectrum fit yields a much lower low energy cut-off, 
around 10 keV, not clearly visible above the large thermal 
component there (Figured). Therefore, we cannot con- 
clude that the true F(E) derived from the primary photon 
spectrum has a well defined low-energy cut-off. 

Figure |H1 gives a further example from the flare of 
August 20, 2002 (Kasparova et al, 2005). We have ob- 
tained regularized F(E) using the method of Kontar et al, 
2005 starting from the total observed spectrum 1(e) and 
the albedo-corrected primary spectrum Ip(e). Note that 
a local maximum in F(E) at ~ 40 keV, difficult to ac- 
commodate in conventional thick-target models (Kontar 
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Fig. 8. Inverted mean electron flux of the August 20, 2002 
solar flare for the time interval 08:25:20- 08:25:40 UT. The 
dash line shows the spectra with albedo correction. The 
confidence intervals represent the range of solutions found 
by allowing the incident photon spectrum to range ran- 
domly within the estimated (instrument + shot noise) er- 
rors. 



and Brown, 2005), is no longer required when albedo is 
properly accounted for. 



6. Discussion and conclusions 

We have shown how the Compton backscatter Green's 
function of Magdziarz and Zdziarski (1995) may be cm- 
ployed to deduce primary hard X-ray spectra from ob- 
servations. We have applied this procedure to RHESSI 
data, particularly from a couple of flares with hard photon 
spectra. Without a treatment of albedo, spectral harden- 
ing found in some flares at photon energies of ~ 40 keV 
appears to require a local minimum in the mean fast elec- 
tron distribution (Piana et al 2003, Kasparova et al, 2005). 
These local minima are particularly interesting, since they 
might, if steep, (Kontar and Brown, 2005) be inconsistent 
with the very widely used collision-dominated thick-target 
model for X-ray production (Brown 1971; Lin and Hudson, 
1976). Here we have seen that a complete treatment of 
albedo removes much of the spectral hardening in this 
photon energy range, potentially restoring the viability of 
the collisional thick target. 

The major assumption here is the isotropy of the down- 
ward directed radiation. At the relevant photon energies 
here 10 — 100 keV, the intrinsic bremsstrahlung cross- 
section polar diagram has a characteristic width of about 
30° (see figure in Massone et al, 2004). Since the emitting 
electron angular distribution will be broadened by pitch- 
angle scattering (Leach and Petrosian, 1981; MacKinnon 
and Craig, 1991), the resulting hard X-ray flux might be 
fairly close to downward isotropic. 



Inclusion of the albedo effect reduces the number of 
energetic electrons required for the production of the ob- 
served spectra. The total flux of energetic electrons 

/•OO 

< vn e nV >= < F{E)nV > dE, (11) 

where v is the average electron velocity. A thin-target fit of 
the observed spectrum of the September, 2002 event gives 

< vn e nV >— 6.2 x 10 54 cm~ 2 sec~ 1 , 15% larger than the 
value obtained for the albedo corrected spectrum . This 
effect is more pronounced in the total energy than in to- 
tal numbers, the energy flux < vEn e nV > being 30% 
overestimated when Compton back-scattering is ignored. 
The influence of albedo can be even more substantial, es- 
pecially for a very flat spectra. For example in the August 
20, 2002 flare analyzed by Kasparova et al. (2005), the ap- 
parent low energy cut-off in observed spectra was found 
around 44 keV, a higher value than normally assumed. 

It should be also stressed here that these are rather 
conservative (lower limit) estimates of the albedo correc- 
tion, assuming an isotropic primary X-ray source. In fact, 
if the electrons are strongly downward directed, the back- 
scattered photons could produce a several times larger 
contribution to the observed spectrum. This would have 
a major effect and clearly one must consider the albedo 
to get a realistic idea of the flare electron spectrum and 
energy budget. The capacity to account completely for 
the effects of albedo also restores some optimism over dis- 
cussing, in terms of X-ray spectra, issues such as the low- 
est energies at which electron acceleration operates (see 
Zhang and Huang, 2004), though electron transport close 
to the thermal speed complicates the interpretation of the 
X-ray spectrum there (Galloway et al., 2005). 
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